The three-dimensional Heisenberg spin glass under a weak random anisotropy. 
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We perform a finite size scaling study of the three-dimensional Heisenberg spin glass in the 
presence of weak random anisotropic interactions, up to sizes L = 32. Anisotropics have a major 
impact on the phase transition. The chiral-glass susceptibility does not diverge, due to a large 
anomalous dimension. It follows that the anisotropic spin-glass belongs to a Universality Class 
different from the isotropic model, which questions the applicability of the chirality scenario. 
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I. INTRODUCTION 

Spin glasses (SG's) are disordered magnetie alloys, 
widely regarded as paradigmatic complex systemsii The 
degree of anisotropy in the magnetic interactions deter- 
mines whether a particular alloy is classified as a Heisen- 
berg or an Ising SG (Ising corresponds to a limit of strong 
anisotropy). Experimentally, anisotropies affect signif- 
icantly the glassy response to external magnetic fields 
and the behavior under cooling protocols.— 

Theorists have privileged the study of the Ising limit, 
in spite of the fact that canonical SG's, e.g. CuMn or 
AgMn, should be rather regarded as Heisenberg, with 
weak anisotropic interactions. Indeed, complications 
arise in the Heisenberg case. In addition to the stan- 
dard SG ordering, Heisenberg systems show as well a 
chiral-glass (CG) phase, where chiralities order— (chirali- 
ties, also named vorticities, reflect the handedness of the 
non-coUincar spin ordering pattern, see definitions be- 
low). 

Probably motivated by failures in early numerical at- 
tempts^ to find a standard SG phase for Heisenberg 
systems, Kawamura proposed a chirality scenario, ex- 
pected to hold for most experimental systemsi^ In the 
ideal, fully isotropic case, the standard SG critical tem- 
perature TsG would be strictly zero, while chiralities 
would order at Tcg > (spin-chirality decoupling). Yet, 
anisotropic interactions (either dipolar, pseudo-dipolar 
or Dzyaloshinskii-Moriyai^) , albeit small, are unavoid- 
able in experimental samples. Hence, the scenario in- 
cludes a decoupling-recoupling hypothesis: weak random 
anisotropic interactions would recouple spins and chiral- 
ities so that Tcg = ^sg > 0. Indeed, the numerical work 
available at the time indicated that very small amounts 
of anisotropy lead to Tgc > Oj^ 

CG ordering may be experimentally investigated 
through the anomalous Hall effect. Due to spin-orbit 
interaction and the spin polarization of the conduction 
electrons, the anomalous Hall resistivity picks contribu- 
tions proportional to the CG order parameter and to its 
corresponding non-linear susceptibilityi^iiS The effective- 
ness of this tool to study non-coplanar orderings has been 



demonstrated in manganites^ii and in a geometrically 
frustrated pyrochlore ferromagnetji^ 

The effect of anisotropies on the critical behavior was 
considered by Bray and Moore^ before the question of 
chiral ordering was raised. They predicted that these 
systems belong to the Ising SG's Universality Class, no 
matter the kind of anisotropic interactions. However, in 
their analysis the assumption was made that Tgc = in 
the isotropic limit (this assumption seemed plausible at 
the time, although we now know that it is incorrect). 

Recent theoretical work has shown that the chirality 
scenario needs some revision. New simulation algorithms 
(allowing to thermalize at lower temperatures than pio- 
neering work^), combined with modern finite-size scal- 
ing (FSS) methods,—"— have provided conclusive evi- 
dence for a standard SG ordering with Tsg > for 
purely isotropic interactions.—"— Only some controversy 
remains on whether Tsg is slightly smaller than Tcg,~ 
or rather the two are compatible within errorsi^ Inter- 
estingly enough, a modern-styled study seems to be still 
lacking for the more realistic case of a Heisenberg SG 
with small random anisotropy. 

Here we show that small anisotropic interactions cause 
that, at variance with the ideal case, the CG suscep- 
tibility no longer diverges at Tcg (i-e. the anomalous 
dimension becomes rjco > 2). In the Renormalization 
Group framework,— anisotropy is a relevant perturba- 
tion. Even if in an experimental sample anisotropies are 
fairly small, the isotropic model is appropriate only for 
moderate correlation length. Closer to the critical tem- 
perature, a new fixed point rules (presumably in the Ising 
SG Universality Class, due to spin- reversal symmetry). 
A slow crossover— from the Heisenberg to the anisotropic 
fixed-point arises upon approaching the phase transition. 
We conjecture that this crossover explains^ experimen- 
tal claims of a non-trivial dependency of critical expo- 
nents on the anisotropy strengthi^i^i Our results follow 
from a FSS analysis of equilibrium Monte Carlo simula- 
tions on system-sizes up to L = 32. Data suggest that 
anisotropies cause a temperature range in which chirali- 
ties order while spins do not (i.e. Tgc < Tcg)- However, 
due to the slow crossover, further research will be needed 
to dismiss spin-chirality recoupling. 
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The remaining part of this work is organized as fol- 
lows. We define the model and describe our numerical 
methods in Sect. |lll We address thermal equilibration, 
a major issue in any spin-glass simulation, in Sect. IIIII 
Our physical results are reported in Sect.|lVl Finally, we 
give our results in Sect. |Vl 



II. MODEL AND SIMULATIONS 
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Since the main types of anisotropic interactions lead to 
the same effective replica Hamiltonian^ it is numerically 
convenient to study short range (pseudo-dipolar) interac- 
tions. Take the Edwards- Anderson model on a cubic lat- 
tice of size L, with periodic boundary conditions. Heisen- 
berg spins occupy the lattice nodes x [S^ = (S*^, S'^, S^), 
Sx ■ Sx = !]■ The Hamiltonian is^ 



H = - 



E 



Sen, 



(1) 



{{x,y): lattice nearest-neighbors). The random 
exchange-couplings, Jxy, arc Gaussian distributed with 

The random Dxy are 3x3 



J. 



x,y 



0, and J'^ y ~ 1. 



DxySy DxySx ' Sy^ 



symmetric matrices (i.e. 5a; 
Their matrix elements are independent and uniformly 
distributed in {^D,D). In most of the work reported 
here D ~ 0.05 (which corresponds to the best studied 
case^), but we will be presenting results for _D = 0.1 as 
weU. 

The ideal limit of a fully isotropic Heisenberg model is 
recovered from Eq. ([T]) by setting D = 0. Once D > 0, 
the original 0(3) symmetry, corresponding to a global 
spin rotation (or reflection), is lost. The only remaining 
symmetry for D > is global spin inversion. 

An instance of the couplings, {Jx,y, Dx%} is named a 
sample. For any physical quantity, we first obtain the 
thermal average, denoted as (...). Only afterwards we 
perform the sample average (denoted by an over line). 

Defining the SG and CG susceptibilities requires real 
replicas. We consider pairs of spin configurations, 5° 
and S*^ , that evolve with independent thermal noise, un- 
der the same couplings and at the same temperature. 
The spin-overlap field is qx = ■ S^., while its Fourier 
transform at wave vector fc, is qsaik) = e'^^'^/iV. 
On the other hand, the local chirality is defined as: 



Cxii Sx 



{Sx X Sx-e^ ) J M — I7 2, 3 , 



(2) 



where is the unit lattice vector along the fi axis. From 
([2]), the chiral overlap-field is k^,^ ~ Cx,^ Cx.^n where 
the superindices a and b correspond to the replicas. Its 
Fourier transform is q'^Q^k) = J2x ^^'^'^ /N. 
The wave-vector dependent susceptibilities are: 



XsG(fc) = N {\qsG{kW) , X'coik) = N (ICgWP)- (3) 



TABLE I: Details of simulations with D = 0.05. For each 
lattice size and temperature, we give the number of simu- 
lated samples. The last row indicates the number of Elemen- 
tary Monte Carlo Steps (EMCS). The L-dependent EMCS 
consisted of 1 Heat-Bath full lattice sweep, followed by 5L/4 
sequential (microcanonical) overrelaxation sweeps. We took 
6 X 10'' measurements per sample, but for L = 32 (15 x 10* 
measurements) . 
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TABLE II: As in Table HI for our simulations with D = 0.1. 



The correlation length, either SG or CG, iai^i^^ 

X(0) ^ 



1 



2sin(fcinin/2) VX(^min) 



1 



(4) 



where fc,nin = (27r/L,0, 0) or permutations 

Our simulation algorithm combines heat-bath with 
microcanonical overrelaxation<21 Both moves general- 
ize straightforwardly to the anisotropic case<^ The 
mixed algorithm is effective for the isotropic Heisen- 
berg s 011 -20,29 g^jj(j fQ]. other frustrated modelsi^ Be- 
sides, we extrapolate to nearby temperatures using a 
bias-corrected^ data reweighting methodi^ Most of our 
simulations were carried out with D = 0.05, see Table [J 
Nevertheless, we did as well some work for D = 0.1, see 
Table HH 



III. EQUILIBRATION 

We considered three thermalization tests. First, con- 
sider the identity (valid for Gaussian-distributed Jx.y)' 



qi 



T 



+ -c/ = o, 



(5) 



where U = - ^i^x,y) Jx,y{Sx ■ Sy)/L^, the link-overlap 
is qi = '2Eix,y){k ■ S'^){Si ■ S'y)/{zL% while q, = 
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FIG. 1: (color online) Sample-averaged A defined in the l.h.s. 
of Eq. ((5| vs. Monte Carlo time, as computed for L = 32 
at T = 0.187 and D = 0.05. The EMCS was defined in 
the caption to Table [H Each point is an average over 3000 
consecutive measurements. 



2J2{x,y}{iSa: ■ Syf)/{zL°) (z = 6 is the lattice coordi- 
nation number). Now, both U and qs equilibrate easily. 
Yet, since qi involves two replicas, it slowly grows from 
zero until its equilibrium value. Thus, a thermalization 
bias shows up as Zi > 0, 17^20,33 r^j^^ time evolution of Z\, 
for L = 32 at the lowest T, is in Fig.[T] Second, we carried 
out the standard logarithmic data binning: we compare 
averages over the second half of the Monte Carlo history, 
with the second fourth, the second eight, and so forth, 
finding stability for three bins. Third, we checked for 
compatibility among reweighting extrapolations for con- 
tiguous temperatures (our simulations at different T are 
statistically independent, see Fig. [2]). 



IV. RESULTS 

Our FSS analysis compares the correlation length in 
units of the lattice size for pairs of lattices {L,2L)i^— 
Dimensionless quantities, such as ^/L, are functions of 
L^/'^iT — Tc), V being the thermal critical exponent. 
Thus, the two curves intersect at Tc{L,2L), see Fig. [5] 
Tc{L,2L) differs from due to scaling corrections (but 
tends to it for large iM-). Our dimensionless quanti- 
ties ^sg/L and ^cg/L, produce two L-dependent critical 
temperatures rsG(-^, 2L) and Tcq{L,2L). We compute 
the anomalous dimensions 77 from the scaling of the sus- 
ceptibilities X [take A; = in Eq. ([3])]. For large L, and 
rj < 2, X diverges as x oc \T — Tc\^ For finite L, 

we consider the susceptibility ratio for xcG and xsG (the 
dots stand for scaling corrections): 



FIG. 2: (color online) Correlation length in units of the lattice 
size vs. T, for all our system sizes a,t D — 0.05. Results for 
both the CG (top), and the SG sectors (bottom). Data 
patches correspond each to an independent simulation (we 
used data reweightingSii) . Inset: CG intersections for pairs 
of sizes (L, 2L). The range of T and (,cg/L differ from main 
plot. 
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TABLE III: Size-dependent critical temperatures Tsg{L,2L) 
and Tcg(L, 2L), and anomalous dimensions 2 — t;sg and 2 — 
?7CG, Eq. ([6]), for the simulations with D — 0.05. Errors were 
obtained with jackknife. 



We discuss first the CG sector. The inset of Fig. [5] 
shows an unusual feature: ^cg/L at the crossing point 
Tcg{L,2L) approaches zero for large L. This is to be 
expected only if 77 > 2>i^ if the susceptibility does not 
diverge at Tc, the correlation length in Eq. ^ scales as 
^/L ~ i-{n-2)/2_ Nevertheless, we still find crossings 
when comparing lattices sizes L and 2L, see Fig. [5] and 
also Ref. [SJ. Crossings are due to the fact that, in the 
large-L limit, the correlation length in Eq. ^ is diver- 
gent in the low-temperature phase. For T < Tc, ^/L 
grows as L^/^ (i.e. the correlation function at large dis- 
tances r goes to a constant with corrections of order 
see e.g. Ref. HI). Yet, the susceptibility ratio in Eq. (|6]) 
is constant for large L, even if 77 > 2. So, 7;cg in Table Hill 
approaches 2 as i grows. 

Besides, it is note worthy that, in spite of the small- 
ness of D, Tcq{L,2L) for D = 0.05 is about twice its 
value for the isotropic model, Tcg{D = 0) « 0.13^-- In 
fact, extrapolating the data in Table Hill as Tcq[L, 2L) = 
TcG + A/L yields Tec « 0.26. 

To further investigate the lacking divergence of xcG at 
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TcG, we consider the integrals^^ 

L/2 

/,=^r^Cp,p(r), (7) 

r-O 

where Cp.p(r) is the plane-to-plane correlation 
functioni^ Note that xcG ~ 2/o, which means 
that plane-to-plane correlation functions decays with r 
slower than the standard point-to-point correlations by a 
factor r^^^. The scaling behavior of the integrals (O is: 
Ik ^ constant in the paramagnetic phase, Ik ^ L^+'^^'^i 
at TcG (if A: + 2 — ?7 > 0, otherwise it is Ik constant), 
and Ik ~ L^'^*' in the CG phase. We show in Fig. [3]-top 
our data for xcG (which is basically 2/o) and, in Fig.[3]- 
bottom, Ii. Note that for T < 0.22 the two integrals are 
diverging with L. On the other hand, for T = 0.22, 0.23, 
Ii grows with L, while XCG does not, as expected for 
2 < ?7cG < 3. 

The behavior of the SG sector is more conventional. A 
remarkable feature in Fig. [2] and Table IIIII is the strong 
scaling corrections in Tsg{L, 2L). We do not consider it 
safe to extrapolate Tsg to its large-L limit, as we are far 
from the asymptotic regime. The SG anomalous dimen- 
sion takes a negative value as L grows (also found in the 
Ising SG, see e.gJ^). 

An intriguing feature is that TgG seems smaller than 
TcG- Indeed, see Fig. d at T « 0.187, where isc/L 
becomes L-independent, S,cg/L is growing fast with L. 
Yet, three caveats prevent us from considering this con- 
clusion as definitive: (i) our lattice sizes are still in a 
strong cross-over regime, hence the final picture could 
change as L grows, (ii) the behavior is rather marginal, 
meaning a larger number of samples would be needed to 
accurately locate Tsg (this is hardly surprising, given the 
large value of exponent i/, and the small 6 exponent, for 
d = 3 Ising SG's), and (iii) when considering a larger 
anisotropy, see below, the effect seems smaller. 

Indeed, we have performed further simulations with 
D = 0.1, up to L = 16. As shown in Fig.|31 the difference 
between Tsg and Tcg is less clearly defined than for D = 
0.05. On the other hand, the chiral-glass susceptibility is 
not divergent at the critical point, in agreement with our 
results for D = 0.05. Consistently with that, the crossing 
points for ^cG /L shift to a smaller height when L grows. 
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FIG. 3: Chiral glass susceptibility xcG (top) and 7i integral 
defined in Eq. [7] (bottom) vs. temperature, for all our system 
sizes with D = 0.05. Lines are guides to eyes. Inset: zoom 
of Ccc/i vs. T 
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T 

FIG. 4: For all our lattice sizes at _D = 0.1, we show the 
chiral-glass susceptibility (top), as well as the CG (center) 
and SG (bottom) correlation lengths in units of the system 
size, as a function of temperature. 



V. CONCLUSIONS 



In summary, we have performed a finite size scaling 
study of the 3d Heisenberg spin glass in the presence 
of a weak random anisotropy, for lattices of size up to 
L = 32. Anisotropics cause that the CG susceptibil- 
ity no longer diverges at Tcg, the chiralities ordering 
temperature. Hence, the anisotropic system belongs to 
a Universality Class different from the isotropic model 
(probably that of Ising SG's). Besides, we found that 
the spin-glass ordering sets up only at Tsg < Tcg- The 
most economic scenario is that actually Tsg = Tcg (the 



apparent difference would be due to finite-size effects). 
In this scenario, chiralities would merely be a compos- 
ite operator (such as, say, the ninth power of the spin 
overlap). However, the would-be intermediate tempera- 
ture region where only chiralities order should be exper- 
imentally detectable through the anomalous Hall effect. 
Numerical studies covering a wider range of values for 
the anisotropic coupling could also help to elucidate the 
situation. 
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